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We have studied far infrared transmission spectra of a'-NaV20s between 3 and 200 cm -1 in 
polarizations of incident light parallel to a, b, and c crystallographic axes in magnetic fields up to 
33 T. The temperature dependence of the transmission spectra was studied close to and below the 
phase transition temperature T c — 34 K. The triplet origin of an excitation at 65.4 cm -1 (8.13 meV) 
is revealed by splitting in the magnetic field. The g-factors for the triplet state are g a = 1.96 ± 0.02, 
gb = 1.975 ± 0.004 and g c — 1.90 ± 0.03. The magnitude of the spin gap at low temperatures is 
found to be magnetic field independent at least up to 33 T. All other infrared-active transitions 
appearing below T c are ascribed to zone-folded phonons. Two different dynamic Dzyaloshinskii- 
Moriya (DM) mechanisms have been discovered that contribute to the oscillator strength of the 
otherwise forbidden singlet to triplet transition. First, the strongest singlet to triplet transition is 
an electric dipole transition where the polarization of the incident light's electric field is parallel 
to the ladder rungs (Ei || a). This electric dipole active transition is allowed by the dynamic 
DM interaction created by a high frequency optical a-axis phonon. Second, in the incident light 
polarization perpendicular to the ladder planes (Ei || c) an enhancement of the singlet to triplet 
transition is observed when the applied magnetic field shifts the singlet to triplet resonance frequency 
to match the 68 cm -1 c-axis phonon energy. The origin of the second mechanism is the dynamic 
DM interaction created by the 68 cm -1 c-axis optical phonon. The strength of the dynamic DM is 
calculated for both mechanisms using the presented theory. 

PACS numbers: 78.30.Hv, 75.10.Pq, 71.70.Gm, 76.30.Fc 



I. INTRODUCTION 

The opening of a spin gap is of fundamental inter- 
est in one-dimensional spin one- half systems. In one- 
dimensional Heisenberg spin chains the coupling between 
the spins and lattice leads to the spin-Peierls instability; 
the atomic distances change together with the nearest 
neighbor exchange coupling between the spins. As a re- 
sult the spin gap opens separating the singlet ground 
state from the excited triplet state. The spin-Peierls 
instability was discovered in organic compounds and 
later on in inorganic CuGeOa 1 . Although a'-NaV205 is 
another spin one-half quasi-one-dimensional compound 
where the spin gap opening and the lattice distortion 
take place simultaneously^, it is different from canoni- 
cal spin-Peierls systems. The magnetic field dependence 
of the phase transition temperature^ T c and the en- 
tropy change^ at T c are not consistent with the magneto- 
elastically driven phase transition in cv'-NaV205, where 
in addition to the displacement of atoms at T c = 34 K a 
new charge order appearsiSuLiii. An extra degree of free- 
dom comes from one electron being shared by the V-O-V 
rung as a'-NaV^Os is a quarter-filled two leg spin ladder 
compound. In the low temperature phase unpaired elec- 



trons on V-O-V rungs shift from the middle of the rung 
to off-center positions and a zigzag pattern of V + and 
V 5+ ions along the ladder legs exists on all ladders2*2*ii. 
There are four ladder planes in the unit cell at ambient 
pressure and under high pressure a series of successive 
phase transitions to phases with more than four planes 
in the unit cell are observed^. 

a'-NaV205 is not an ideal one-dimensional spin system 
as revealed by inelastic neutron scattering (INS)iii4. In 
addition to the dispersion of magnetic excitations along 
the ladder direction two dispersion curves with a rather 
small dispersion of 1.2 meV along the ladder rung direc- 
tion are observed. One curve is at 8meV and the other 
at lOmeV at the center of the Brillouin zone. More exact 
spin gap value, 8.13 meV (65.5 cm -1 ), has been deter- 
mined by the high field electron spin resonanc o 15 i 16 . 

Doubling of the lattice constants along a- and b-axes 
and quadrupling along c-axis create additional Raman 
and infrared active modes at the phase transition tem- 
perature. The question is whether they are all zone- 
folded lattice modes or some of them are spin excita- 
tions. The controversial modes are infrared active modes 
polarized along the c-axis at 68 and 106 cm -1 and Ra- 
man modes with frequencies nearly matching the fre- 
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FIG. 1: Two spins coupled by the isotropic exchange coupling, 
J = 1, and the DM interaction, D || Bo, D = 0.4. (a) - energy 
levels; |t_) = |T_) and |t+) = |T+) are pure triplet states in 
any field, (b) - non-zero transition probabilities If from the 
ground state |s) to the triplet state |tj) for a given orientations 
i = x or y of the alternating magnetic field Hi . 



FIG. 2: Two spins coupled by the isotropic exchange coupling, 
J = 1, and the DM interaction, D 1 B ,D = 0.4. (a) - 
energy levels; |to) = |To) is a pure triplet state in any field, 
(b) - non-zero transition probabilities I\ from the ground 
state |s) to the triplet state |tj) for given orientations i = 
x,y,or z of the alternating magnetic field Hi. 



quencies of the two infrared modes. The origin of the 
68 and 106 cm -1 excitations is of fundamental interest 
since the 68 cm -1 excitation is almost degenerate with 
the 65.5 cm -1 triplet. Spin chain models do not predict a 
bound singlet excitation being degenerate with the triplet 
excitation 17 i 18 i 19 . 

Electric and magnetic dipole transitions between sin- 
glet and triplet states are forbidden in principle. The 
reason for this is the different parity of the ground sin- 
glet state and the excited triplet state. The singlet state 
is anti-symmetric and the triplet state is symmetric rel- 
ative to the interchange of two spins. Electric dipole 
or magnetic dipole operators, responsible for the optical 
absorption, will not couple these two states. An anti- 
symmetric interaction can mix singlet and triplet states. 
If such interaction exists the transitions are partially al- 
lowed. The strength of the partially allowed optical tran- 
sition, either magnetic dipole or electric dipole, depends 
on the orientation of the light polarization and applied 
magnetic field with respect to the crystal axes. Using 
infrared spectroscopy it is possible not only to extract a 
separation of energy levels in a spin system, but also find 
the origin of the transition, either electric dipole or mag- 
netic dipole, and the orientation of the anti-symmetric 
interaction. 

In this paper we study spin gap excitations and 
phonons in a'-NaV20s using far-infrared spectroscopy. 
To explain the experimentally observed singlet to triplet 



absorption we present a calculation of the dynamic 
Dzyaloshinskii-Moriya (DM) absorption mechanism by 
numerical diagonalization of the spin-phonon hamilto- 
nian. Analytical results in the perturbation theory for 
the small dynamic DM interaction are given. Also, the 
theory of the second relevant mechanism, the magnetic 
dipole active static DM mechanism is presented. 



II. STATIC AND DYNAMIC 
DZYALOSHINSKII-MORIYA INTERACTION 

The anti-symmetric DM interaction, introduced by 
Dzyaloshinski: 20 and Moriya2I, is a combination of su- 
perexchange and spin-orbital interactions and is linear in 
spin-orbital coupling. For a particular spin system the al- 
lowed components of the DM interaction arc determined 
by the symmetry of the spin comple:»22i2i. Corrections 
to the energy spectrum due to the DM interaction are 
usually small because the correction is proportional to 
D 2 /A, where D is the magnitude of the DM interac- 
tion and A is the separation of the singlet and triplet 
energy levels. In addition to the DM interaction there 
is a symmetric spin-spin interaction that is quadratic in 
spin-orbital coupling. Although the symmetric interac- 
tion does not couple the singlet and triplet states it affects 
the splitting of the triplet state sublevels. It was shown 
by Shekhtman et alm2£& for a single bond superexchange 
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that the triplet state remains degenerate in zero magnetic 
field if both symmetric and anti-symmetric spin-spin in- 
teractions are taken into account. 

Although the corrections to the energy levels are small, 
the mixing of the singlet and the triplet state by the anti- 
symmetric interaction could be enough to produce exper- 
imentally detectable optical singlet to triplet transitions. 
The transition probabilities for the magnetic dipole op- 
erator in spin chains with DM interactions were calcu- 
lated by Sakai et al£*. The idea that the electric dipole 
singlet to triplet transition is partially allowed when an 
optically active phonon lowers the crystal symmetry and 
therefore creates a dynamic DM interaction, was put for- 
ward by Cepas et al£^2&. Below we calculate the en- 
ergy spectra and transition probabilities using a simple 
two spin model. In the case of the static DM mechanism 
the full spin hamiltonian with Shekhtman corrections is 
used. For the dynamic DM interaction we extend the 
theory of Cepas et al. beyond the perturbation theory 
and solve the hamiltonian by exact diagonalization. 



where Hi is the magnetic field component of the light, 
g is the electron g- factor, and /^s is the Bohr magneton. 
The zero-order hamiltonian is 

H = JS t ■ S 2 + ,g /is Bo ■ (Si + S 2 ), (2) 

where J is the isotropic exchange coupling between spins 
Si and S 2 , and B is the applied static magnetic field. 
The first and second order corrections are (Eq.(2.19) from 

H sx = -^S 1 -S 2 + ^jS 1 -r>T>-S 2 + H DM , (3) 

where D is the DM vector and we have separated the 
anti-symmetric DM interaction 

H DM = D- [Si x S 2 ]. (4) 



A. Magnetic dipole transitions and 
Dzyaloshinskii-Moriya interaction for a spin pair 

In this section we calculate the energy spectrum 
and the magnetic dipole transition probabilities for two 
isotropically exchange coupled spins (S — 1/2) in the 
presence of the anti-symmetric DM interaction and the 
second-order symmetric interaction. We find the eigen- 
values and eigenstates of the hamiltonian H — Hq + H sx 
and calculate the transition probabilities from the ground 
state induced by the magnetic dipole operator 

H m d = 5/isHi • (Si + S 2 ), (1) 



We choose singlet and triplet as the basis of eigenstates 
since they are the eigenstates of Hq. These states are 
the singlet |S) = (| + -) - | - +>)/\/2 and the three 

components of the triplet |T_) = | ), |T ) = (| + -) + 

| — h))/-\/2, and |T + ) = | + +). The spin quantization 
axis z is chosen parallel to the applied field Bo and for 
a single spin (+|5 2 |+) = -(-|&|-) = 1/2. Let the 
state vector be |*) = (T + , T , T_, S). We diagonalize 
the hamiltonian H = Hq + H sx for two orientations of 
applied field, along the DM vector and perpendicular to 
the DM vector, denoting the eigenstates by (t+, to, t_,s). 

Bo || D || z. The hamiltonian in the matrix represen- 
tation is 
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where G z = g/is Bq is the Zeeman term. For arbitrarily 
chosen J = 1 and D = 0.4 the energy levels are shown 
in Fig.^a). The symmetric part of H sx adds a correc- 
tion D 2 /(16J) to all energy levels except |T ) where it is 
— 3D 2 /(16J). This correction for |To) is partially can- 
celed by Hdm that mixes |S) and |T ). As a result 
the triplet state sublevels stay degenerate in zero mag- 
netic field as was pointed out in refi2Si2i. The net effect 
of H sx in zero field is to lower the singlet state energy 
by — 3D 2 / (16 J) and to raise the triplet state energy by 
7J 2 /(16J). 



The transition probability for the magnetic dipole op- 
erator (JIJ from the ground state is calculated as If — 
\(tj\Sn + S^s)! 2 , i = x,y,z. The alternating magnetic 
field Hi polarized along x or y axis (perpendicular to Bo 
and D) gives non-zero intensities as shown in Fig.^b). 
This is because |s) has the triplet component |To) mixed 
in and the transitions from |To) to |T_) and |T+) are 
allowed by and S y operators. The transition proba- 
bilities do not depend on the strength of the applied field 
since the mixing of \S) and |T ) is independent of B . 

Bq ID j) y. The hamiltonian is 
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In this field orientation |T_) and |T + ) are mixed into |S) 
by Hum and | to ) remains a pure state, |to) = |To). Note 
that there is an avoided crossing at g^BBo ~ 1 between 
|s) and |t_) as shown in Fig.[2[a). 

The transitions from |s) to |t_) and |t+) are observed 
when Hi || Bo, see I~ and 7+ in Fig.EJb). In high mag- 
netic field I~ prevails over 7+ because the mixing of |T_) 
into the ground state increases and the mixing of |T + ) 
decreases. Finite transition probability 7° to the |to) is 
observed in small fields when Hi _L Bo, D whereas 7° = 
(Hi || D) as Bq approaches zero. Both transition proba- 
bilities are determined by the amount |T_) and |T + ) are 
mixed into the ground state since transition operators 
and S y couple these two states to the |to) = |T ) state. 
7° and I® gain intensity as the ground state changes into 
|T_) with increasing field. 

In summary, the following selection rules are observed 
for the magnetic dipole transition from the singlet to the 
triplet state in the presence of DM interaction. First, if 
the magnetic field is parallel to the DM vector D || Bo, 
transitions to the triplet state sublevels |t_) and |t+) arc 
observed. These transitions have field-independent in- 
tensities and do not depend on polarization in the plane 
perpendicular to the DM vector, Hi _L D . Sec- 
ond, if the magnetic field is perpendicular to the DM 
vector Bo -L D, then in small fields (Bq <C J/g^s) 
the transition to |to) has a weak field dependence and 
is observed in polarization Hi _L Bo,D. The transi- 
tions to |t_) and |t+) are observed in Hi || Bo polariza- 
tion. In this polarization in magnetic fields, g/iBBo > D, 
the transition probability to |t_) increases and to |t + ) 
decreases with increasing field. Sakai et al 2A calcu- 
lated magnetic dipole transition probabilities for inter- 
acting spin chains using a 16 spin cluster. In their model 
the symmetric anisotropic superexchange was not con- 
sidered. Our single dimer model leaves the triplet levels 
degenerate whereas the degeneracy is lifted in their cal- 
culation. Whether the degeneracy will be lifted or not 
when the symmetric anisotropic superexchange in addi- 
tion to the anti-symmetric DM interaction is included in 
their model needs a separate study. 



B. Electric dipole transitions and dynamic 
Dzyaloshinskii-Moriya interaction for a spin pair 

We show that the electric field component of the far- 
infrared light Ei that couples to an optically active 
phonon can cause transitions between singlet and triplet 



states if this phonon creates a DM interaction by lattice 
deformation. 

Electric dipole coupling between the phonon and the 
light in the long wavelength limit is 

V = eQE 1 , (7) 

where e is an effective charge associated with a lattice 
normal coordinate Q. Here we assume that the electric 
field is polarized parallel to the electric dipole moment of 
the normal coordinate Q and we have dropped the time 
dependence of V. 

We expand the DM vector D into a power series of Q 

D(Q) = D(0) + !^| Q=o Q + .... (8) 

The first term is the static DM vector in the absence of 
lattice deformation. We already demonstrated in the pre- 
vious section that this interaction gives rise to magnetic 
dipole transitions between singlet and triplet states. Here 
for simplicity we take D(0) = 0. We will ignore terms 
quadratic in D in H sx (Eq[3J) because these symmetric 
interactions will not give us any transitions between sin- 
glet and triplet states. Leaving out higher order terms of 
Q we get for the DM interaction 

H DMQ = QD Q • [Si x S 2 ], (9) 

where Dq = ^ | g =0 . For the phonon system we use the 
secondary quantization presentation. The lattice normal 
coordinate Q can be presented in terms of creation and 
annihilation operators a' and a, Q — q{a* +a), where q is 
the transformation coefficient. Since we left out Q 2 and 
higher order terms in Eq.(JHJl, the dynamic DM interaction 
will couple two phonon states \n) and \n'}, where n' = 
n± 1; n is the occupation number of phonons in mode Q. 
The hamiltonian of the coupled spin-phonon system is 

H = ^cj p a t a+JSrS 2 +.gAiBBo-(Si+S 2 )+i7{ 3 1 ], /Q , (10) 

77^ Q = g (at + a)D Q -[S 1 xS 2 ], (11) 

huip is the phonon energy. In the low temperature limit 
ksT Tiuj the thermal population of phonon states is 
low, (n) «0, and we can consider only the phonon states 
with either or 1 phonon, |0) and |1). After diagonal- 
ization of the hamiltonian (|l(Jfl we treat V , Eq.Q, as a 
time-dependent perturbation to calculate the transition 
probabilities from the ground state. We choose |S) and 
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FIG. 3: Two spins coupled by the isotropic exchange coupling, 
J — 1, and the dynamic DM interaction, qDq = 0.4, created 
by the phonon with a frequency Tiuj v = 1.2; Dq || Bo. (a) - 
energy levels; |t±0) and |t±l) are pure states |T±0) and |T±1) 
in any field, (b) - non-zero electric dipole transition probabil- 
ities from the ground state to the coupled spin-phonon state 
|toO) and |sl). The graph is not extended above the field 
where the ground state changes from the singlet |s0) to the 
triplet |t_0). 



FIG. 4: Two spins coupled by the isotropic exchange coupling, 
J — 1, and the dynamic DM interaction, Q-Dq = 0.4, created 
by the phonon with a frequency hu) p = 1.2; Dq _L Bq. (a) - 
energy levels; |to0) and |tol) are pure states |To0) and |Tol) in 
any field, (b) - non-zero electric dipole transition probabilities 
from the ground state to the coupled spin-phonon state |t_0) 
(dotted), |t+0) (dash-dotted), and |sl) (solid). The graph is 
not extended above the field where the ground state changes 
from the singlet |s0) to the triplet |t_0). 



|Tj) with the quantization axis along the applied field Bo 
as the basis for the spin states. Let the state vector be 
* = (T + l,T + 0,T l,ToO,T_l,T_0,Sl,SO). The hamil- 
tonian is diagonal in this basis except for the last term, 
Hjjmq ■ The eigenstates are labeled by |sn) and |t»n) 
where n = or 1. We solve two separate cases, Dq || Bo 
and Dq _L B . 

Dq I) Bo || z. In this field orientation Dq = 
(0,0, Dq). The diagonal elements are the same as in 
except that the phonon energy fiu> p will be added 
if n = 1. Beside the diagonal elements the non- 
zero elements of the hamiltonian i|lU|) are the ones ere- 

ated by H$ MQ : (Sl|^ o |T 0) = (S0|^ o |T l) = 

-(T l|if^ Q |S0) = -(T 0|if^ Q |Sl) = iqD Q /2. The 
energy levels, calculated for J = 1, Tiuj p = 1.2, and 
qDQ — 0.4, are shown in Fig.|3Ja). The largest repulsion 
is between | si) and |trj0), which are the linear combina- 
tions of |S1) and |To0). The other two mixed together 
states are | SO) and |T 1) giving us the ground state |s0) 
and |tol). One has to keep in mind that not only the spin 
states are mixed, but also the phonon states |0) and |1) 
are mixed. All other 4 states that involve triplet states 
|T±) are pure states. 

The splitting of energy levels has to be taken with some 



precaution. The splitting due to the dynamic DM is ob- 
served when there is one phonon excited, n = 1. This 
is not the case at thermal equilibrium at low T when 
(n) = (the possible role of zero-point vibrations is ig- 
nored in our approach). If the phonon is brought to the 
state n = 1 by the light-phonon interaction J7J the ef- 
fect of one phonon on the shift of energy levels should 
be observed in the experiment, unless it is much smaller 
than the lifetime or inhomogeneous broadening of energy 
levels. The magnitude of the shift and whether it could 
be observed in the experiment or not will not affect con- 
clusions about the transition probabilities. 

Calculation of the transition probability | (tjn'|y|s0) | 2 
is straightforward since V couples states that are diag- 
onal in the basis of pure spin states and non-diagonal 
in the basis of phonon states |0) and |1). Two transi- 
tions from the ground state |s0) are observed, to | si) and 
|toO), shown in Fig.HJb). If the dynamic DM is zero, 
then |s0) = |S0), |sl) = |S1), and |t 0) = |T 0>. The 
transition from |s0) to | si) is an ordinary absorption of 
an infrared photon fiu) p = E s \ — E s q with probability 1, 
and the transition to |toO) has zero probability. When the 
dynamic DM interaction is turned on, additional absorp- 
tion sets in and a photon of energy E to o — E S Q is absorbed. 
This can be viewed as a virtual excitation of a phonon by 
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the light to the state |1) while the spin state remains sin- 
glet, and then the dynamic DM interaction @ brings the 
(virtual) phonon back to |0) while changing the spin state 
to | To). The polarization of the absorbed photon with 
respect to the crystal axes is determined by the phonon 
states involved. 

Dq _L B II z. We take D Q = (Q,D Q ,Q). Be- 
side diagonal elements there are 8 non-zero elements 



T + 1\H^ MQ \S0), (T + 0\H 



(1) ISIS 
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DMQ 
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(T_0|if^ Q |Sl), (Sl|ff^ Q |T + 0), (Sl|F^ o |T_0), 
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DMQ I 1 + 
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T_l), all equal to 



^/2qDq/A. The energy levels, calculated for J = 1, 
huj p = 1.2, and qDq — 0.4, are plotted in Fig.^Ja). 

The strongest mixing occurs between | SI) and |T±0) 
giving the eigenstates | si) and |t±0). Also, there is an ad- 
ditional mixing between |T_0) and |T + 0) levels in small 
fields. 

In small fields, gfisBo < J, the mixing of | SO) and 
|T±1) is less pronounced since their separation is larger 
than the separation of |S1) and |T±0). Therefore, for the 
analysis of the transition probabilities the ground state 
can be taken as pure | SO) and the transition probabilities 
are mainly determined by the mixing between |S1) and 
|T±0). The effect of mixing of |T±1) into the ground 
state has a secondary effect on the transition probabili- 
ties. Transitions from the ground state to three excited 

states |t+0), 1 1 0) , and |sl) have non-zero probabilities 

(see Fig.0Jb)). The transition probability to the | si) 
state increases with field because | si) changes gradually 
from the mixed state into |S1). The transition probabili- 
ties to 1 0) and |t_0) state decrease as the field increases 
because the amount of |S1) mixed into them decreases. 
Again, the polarization of the absorbed photon with re- 
spect to the crystal axes is determined by the phonon 
states involved. 

Perturbation theory. Analytical results can be ob- 
tained in the limit |i?si — E^^l ^ qDq, i — — ,0,+. 
This case holds when \huj p — (J ± gfisBo)] ^S> <Z-Dq- We 
find the first order perturbation corrections to the states 
|Sn) and iT^n), where n = 0, 1, using H^ M q, Eq. (|11|) . as 
perturbation. Then the transition probabilities are cal- 
culated between the new states |s0) and 1 0) as was done 
in the exact treatment. 

For Dq II B || z we get 



(12) 



where I p — (eqEi) 2 is the light absorption intensity by 
the infrared-active phonon. The transition probability 
from | So) to the triplet level | to ) is independent of the 
magnetic field. 

For the perpendicular case, Dq _L Bo || z 

(qD Q f{hu p f 



|<t_0j^|s0)|^ = 7 3 



|(t + 0|V|s0)| 2 



p 2[{hLo p y-{j-gLL B B yf 

(qD Q ) 2 (hw p ) 2 
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(14) 



If huj p <C J then the intensity of the transition from |s0) 

to 1 1 0) increases with the magnetic field and decreases 

for the transition to |t + 0). If huj p 3> J then the intensity 
of the |s0) to |t+0) transition increases and of |s0) to 

1 1 0) decreases. In the perturbation limit the zero field 

intensities of the transitions from |s0) to 1 1 0) and |t+0) 

are equal. 

In summary, the following selection rules are obtained 
for the electric dipole transition from the singlet to the 
triplet state in the presence of the dynamic DM interac- 
tion. 1. The polarization of the transition: E x is parallel 
to the dipole moment of the optically active phonon that 
creates the dynamic DM interaction. 2. The orientation 
of the dynamic DM vector Dq is determined by the sym- 
metry of the lattice distortion created by the optically 
active phonon. 3. If Bo || Dq a magnetic field indepen- 
dent transition probability to the triplet state sublevel 
|to0) is observed. 4. If Bo J- Dq magnetic field depen- 
dent transition probabilities to the triplet state sublevels 
|t + 0) and 1 1_ 0) are observed. 



III. EXPERIMENTAL 

We studied several single crystals of a'-Na. x V20 y from 
the batch El 0621. According to the heat capacity mea- 
surements these crystals have T c — 33.9 K, and the chem- 
ical composition x = 1.02 and y = 5.06. The (ab)-plane 
properties, Ei_Lc, k || c, were studied on three single 
crystals, 600 /im (area in the (ab)-plane 21 mm 2 ), 120 /im 
(20 mm 2 ) and 40 (im (3.5 mm 2 ) thick in c-direction. The 
(bc)-plane properties, Ei_La, k || a, were measured on 
a mosaic of three crystals, each approximately 650 ^m 
thick, with a total area of 19 mm 2 in the (bc)-plane. The 
(ac)-plane properties, Ei_Lb, k || b, were measured on 
a mosaic of 7 crystals, each approximately 800 /im thick, 
with a total area of 8.5 mm 2 in the (ac)-plane. 

The far-infrared measurements were done with a polar- 
izing Martin-Puplett Fourier transform spectromete»2£. 
Light pipes were used to guide the far infrared light into 
the sample cryostat equipped with a 12 T Oxford Instru- 
ments magnet and two silicon bolometers from Infrared 
Laboratories operated at 0.3 K. A rotatable polarizer was 
mounted at the end of the light pipe in front of the sam- 
ple to control the polarization of light. Spectra were 
recorded at 0.2 to 0.3 cm -1 resolution. The magnetic 
field was applied parallel to the direction of light propa- 
gation (Faraday configuration, k || Bo) or perpendicular 
to the light propagation (Voigt configuration, k_LBo). 
Measurements above 12 T were performed at the Na- 
tional High Magnetic Field Laboratory on a 33 T Bit- 
ter magnet in the Faraday configuration using a Bruker 
IFS 113v infrared spectrometer and a 4K silicon bolome- 
ter from Infrared Laboratories. 

The anisotropic power absorption coefficient ai{ui) 
(i = a, b, c) was calculated from the measured trans- 
mission Ti{oj) assuming one back reflection from the 
crystal front face and one from the back face, Ti(uj) = 



7 



TABLE I: Absorption line positions wo (cm -1 ), full widths at 
half maximum 7(cm _1 ), and areas (cm -2 ) in 4.4 K and 40 K 
spectra of a'-NaV2 0s in zero magnetic field. Index T refers to 
the singlet to triplet excitation and F to the Fano line shape. 
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FIG. 5: Absorption spectra of a'-NaV 2 5 below, T = 4.4 K, 
and above, T = 40 K, the phase transition temperature, (a) 
- Ei || a, (b) - Ei || b. The inset shows the singlet to triplet 
excitation at 65 cm -1 . 



IV. RESULTS 



(1 — Ri) 2 exp[— ati(u))d\, where d is the thickness of the 
crystal. We used a frequency independent value for the 
reflectance coefficient Ri = [(rii — l)/{rii + l)] 2 . In- 
dexes of refraction, rii, at 4K and teraherz frequencies 
are n a = 3.64, n b = 3.16 and n c = 2.70 (Ref3). Accord- 
ing to another paper 30 indexes of refraction at 4K and 
0.55 cm- 1 are n a = 3.07, n b = 3.19 and n c = 2.03. We 
calculated n a jn b = 1.16 from the fringe pattern in the 10 
to 50 cm -1 range using our transmission data. This ratio 
is more similar to the ratio n a /n b — 1.15 from RefiS and 
therefore we used their values for refraction indexes. The 
real part of the conductivity in units of f2 _1 cm -1 is 

ffi(w)=T*(a;)ai(a;)/(12(]7r), (15) 

where in the limit of weak absorption we take rii to be 
independent of frequency. 



A. Absorption spectra and their temperature 
dependence 

The absorption spectra at temperatures above and be- 
low the phase transition temperature T c for incident light 
polarizations Ei || a and Ei || b are shown in Fig.JSJand for 
Ei || c in Fig. |SJ Below T c several new lines appear. The 
line parameters are listed in Table |l| The full width at 
half maximum (FWHM) of the narrowest lines is deter- 
mined by the used instrument resolution, 0.2 cm -1 . The 
best fit for the resolution limited narrow lines was ob- 
tained using gaussian line shapes, otherwise a lorentzian 
line shape was used. 

The Ei || a absorption at 40 K is dominated by a con- 
tinuous absorption steadily increasing towards high fre- 
quencies. Above 180 cm -1 the absorption is too strong 
and our data is not reliable above this frequency. There 
are two derivative-like absorption lines, one at 91.2 cm -1 
and the other at 140 cm -1 , the latter being relatively 
broad and has its phase opposite to the 91.2 cm -1 line 
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FIG. 6: Absorption spectra at 4.4 K and 40 K in Ei || c polar- 
ization. 
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FIG. 7: Temperature dependence of the singlet to triplet ab- 
sorption spectrum in Ei || a polarization. 
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FIG. 8: Zero magnetic field temperature dependence of the 
normalized absorption line area (panel a), normalized reso- 
nance frequency (b), and FWHM (c) for the singlet to triplet 
transition at 65.4 cm -1 Ei || a (circles), and for the 68 cm -1 
c-axis phonon (crosses). Inset to (b): the solid line is a fit of 
the S to T transition energy A(T)/A(4.7K) = (1 - T /T c )> 3 
above 20 K; T c = 33.9 K, (3 = 0.039 ± 0.002. Additionally the 
X-ray diffraction peak intensity from Ref»S is plotted with a 
dashed line in panel (a) and the normalized gap at 9.8 meV 
measured by INS- is shown with triangles in (b). 



phase. There is a narrow line at 137.6 cm" 1 on top of 
the 140 cm" 1 line. When T is lowered, the absorption 
continuum is diminished over the entire frequency range. 
Still, substantial absorption continuum remains above 
130cm- 1 . At T = 4K the 91.2 and 140cm" 1 lines have 
a normal absorption-like line shape. New lines appear 
as well: doublets at 102, 127, 147 cm" 1 and single lines 
at 112, 157, 168cm -1 . The inset in Fig. 03a) shows a 
weak line at 65.4 cm -1 , the singlet to triplet excitation. 
The magnetic properties of this transition are described 
in more detail in section llVBI The T dependent spectra 
of the singlet to triplet transition are plotted in Fig.0 
To better extract this relatively weak line the difference 
of two spectra at fixed T, one measured in T field and 
the other measured in 5 or 10 T field, was calculated. 
One can see that the negative line, representing the zero 
field spectrum, together with the positive features (the 
lines in 5 or 10 T spectra) shift to lower frequencies as 



the temperature increases and at the same time the lines 
lose intensity and broaden. The temperature dependence 
of line parameters is plotted in Fig.|5] 

The Ei || b absorption at 40 K is dominated by a strong 
phonon line at 180 cm" 1 and a weaker line at 215 cm" 1 , 
Fig. Efb) . The temperature-dependent absorption con- 
tinuum, as was observed in Ei || a polarization, is absent 
in Ei || b (note the different vertical scales in Fig. [3^ and 
Eb). Several lines, at 101.4, 111.7, 126.7, and 127.5cm" 1 , 
have the same frequency as in the a-axis spectrum, al- 
though different strength. The line parameters are listed 
in Tabic [I] The temperature evolution of the multiplct 
of seven lines in the Ei || b spectrum below 40 cm" 1 is 
shown in more detail in Fig.® Below T c = 34 K a broad 
line appears at 32 cm -1 . Another line appears at 26 cm -1 
as T is lowered. At low T those two features split into 
a doublet and a quintet. The line at the highest fre- 
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FIG. 9: Temperature dependence of seven low frequency b- 
axis phonons. 
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FIG. 10: The temperature dependence of the 68 cm -1 c-axis 
phonon absorption spectrum. 



quency, 39 cm -1 , is relatively weak compared to other 
six lines. The spectra shown in Fig.^Hhave not been cor- 
rected for light interference fringes in the sample. The 
60 K spectrum could be used as a background, but with 
some precaution as some intensity is lost between 40 and 
70 cm -1 when T is lowered. 

The Ei || c absorption spectrum has no sharp fea- 
tures below 140 cm" 1 in the high temperature phase at 
T = 40 K (Fig. El). When T is lowered below 34 K new 
lines evolve at 68, 106 cm -1 , and two doublets around 
124 and 132 cm -1 . Since the low T transmission was too 
small at the transmission minimums and the absorption 
coefficient cannot be determined accurately, only the up- 
per limits for the line widths and the lower limits for the 
line areas are given in the TableQ] The temperature de- 
pendence of the 68 cm -1 line is shown in Fig.^J The 
line shifts to lower frequency and broadens as T is raised 
from 4.5 K. Some intensity change is still observed above 
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FIG. 11: Magnetic field dependence of the singlet to triplet 
transition spectra in Ei || a polarization. Spectra have been 
shifted in vertical direction. Measurements were done in Voigt 
(a, b) and in Faraday configuration (c, d). 



T c between 34 and 35 K, but there are no visible differ- 
ences between the 35 K and 37 K spectra. The temper- 
ature dependence of line parameters is plotted in Fig-IHl 
together with the singlet to triplet resonance data. Since 
the largest experimentally detectable absorption was lim- 
ited by the thickness of available crystals the line area has 
been reliably determined only above 25 K. Below 25 K the 
plotted FWHM is the upper limit for the line width. 
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FIG. 12: The magnetic field, Bo || b, dependence of the singlet 
to triplet resonance line area at 4.4 K in Ei || a polarization. 
The values for the transition to ms = — 1 are plotted in the 
negative field direction and for the transition to ms = 1 in the 
positive field direction; g a ^B = 0.922 cm -1 /T. The lines show 
the theoretical transition probability for the electric dipole 
transition for a set of parameters where u p is the resonance 
frequency of an optical phonon coupled to the spin system 
by the dynamic DM interaction: solid line - lo p — 180 cm -1 , 
dashed line - lu v = 518 cm -1 . 



B. Singlet to triplet absorption in magnetic field 

We studied the magnetic field effect on the absorption 
spectra below 130 cm -1 at 4.4 K in magnetic fields up to 
12 T in all three polarizations, Ei |j a, Ei || b, and Ei || c. 
We did not see any line shifts nor intensity changes except 
for the 65.4 cm -1 absorption line. 

For Ei || a polarization the magnetic field effect on the 
65.4 cm -1 line is demonstrated in Fig. ^2 When the field 
Bo is parallel to the a- or b-axis, the line splits (panels 
(a), (b), and (c)). In Bo || c configuration the line does not 
split or change its intensity (panel (d)). From the field 
dependence of the resonance frequency the 65.4 cm -1 ab- 
sorption line can be identified as a transition from a sin- 
glet ground state, S = 0, to a triplet excited state, S = 1 
Light is absorbed, depending on its polarization, either 
by transitions to ms — 1 and ms = — 1 triplet levels 
when Bo J- c or to the magnetic field independent ms = 
level when Bo || c. The measurements were extended to 
33 T in one polarization and field orientation, Ei || a, 
Bo || b. The field dependence of the singlet to triplet 
transition line areas is shown in Fig.^] The lines were 
fitted with a lorentzian function, full width at half max- 
imum (FWHM) being between 0.5 and 0.6 cm -1 . The 
line areas of the transitions to the ms = — 1 and 7715 = 1 
levels depend weakly on the magnetic field. Below (see 
section IY A 1(1 we calculate the electric dipole transition 
probability using the presented theory of the dynamic 
DM effect, and compare it to our measurement results. 
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FIG. 13: Differential absorption spectra of the singlet to 
triplet transition at 4.4 K in Ei || a and Ei || c polarizations. 
The Ei || c spectrum has been multiplied by 10 and offset in 
the vertical direction. Both spectra have been measured in 
Bo || a, Hi || b geometry. 

An important question to answer is which component 
of light, Ei or Hi, interacts with the spin system. What 
is common for the data presented in Fig.^2 is that all 
these measurements were done with the light polarized 
along the ladder rungs, Ei |j a. We made complementary 
measurements rotating the incident light's polarization 
by 90°, thus interchanging the orientations of Ei and 
Hi, and found that the singlet to triplet transition was 
at least ten times weaker if Ei _L a. While the data 
presented in Fig.lTTl panels (a), (b), and (d) still leaves 
open the possibility that the singlet to triplet transition 
could be a magnetic dipole transition where Hi || b, there 
are two other observations that contradict this. First, in 
Fig-El ( c ) Hi is not parallel to the b-axis while Ei is still 
parallel to the a-axis and the transition is still strong. 
Second, in Fig. 1131 two spectra with Hi || b are presented 
for the Bo || a applied magnetic field orientation. The 
upper curve with Ei || c shows about 10 times weaker 
absorption on the singlet to triplet transition than the 
lower curve with Ei || a. Since the same orientation of 
Hi gives different intensities the orientation of magnetic 
field component of light is irrelevant. Summarizing the 
results presented in Fig.llllandll3lwe conclude that the 
strongest contribution to the singlet to triplet absorption 
comes from an electric dipole transition with the dipole 
moment along the a-axis. 

There are other contributions to the singlet to triplet 
absorption, although much weaker than in Ei || a polar- 
ization, as is shown in Fig. 1131 The upper spectrum is 
for the c-axis polarized light. Two vertical dashed lines 
mark the area where the strong absorption due to the 
68 cm -1 excitation takes place in this polarization (see 
also Fig. E| . One can see that the absorption lines closer 
to the 68 cm -1 excitation are stronger than the lines fur- 
ther away. 
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The line areas measured in different geometries as a 
function of magnetic field strength are plotted in Fig.1141 
We used the method of differential absorption, where 
spectra taken in different magnetic fields are subtracted 
from each other to detect weak transitions. In this differ- 
ential method the transition from the singlet state to the 
ms — triplet state escapes detection (unless the inten- 
sity depends on the field) since the energy of this triplet 
level does not change with magnetic field. Therefore, 
only the intensities of the transitions to the rug = — 1 
and nis = 1 triplet levels can be detected and are plot- 
ted in Fig. 1141 An interesting finding is the enhancement 
of the singlet to triplet transition in Ei |j c polarization 
close to the 68 cm -1 excitation. It is natural to associate 
the oscillator strength of this weak transition with the in- 
teraction between the spins and the 68 cm -1 excitation, 
which has a dipole moment along the c-axis. The excep- 
tion is Bo || b orientation where there is no enhancement 
in Ei || c polarization. In this field direction only the 
transition to the ms = level has a non-zero oscillator 
strength but is not detected because of the measurement 
method. 

Besides the strong Ei || a absorption and the res- 
onantly enhanced Ei || c absorption there is a field- 
independent oscillator strength, as is seen in Fig. lMf a). 
No polarization anisotropy is observed there in contrary 
to the first two mechanisms of singlet to triplet transi- 
tions that can be recognized by their polarization depen- 
dence. In high fields the transitions to the ms = — 1 
and ms = 1 levels have the same strength in Ei || b 
and Ei || c polarizations (Fig. ll4f a)). When Bo || c and 
Ei || b fFig. lMf b'O the intensity of the transitions to the 
ms = — 1 and ms = 1 levels is zero and in this configura- 
tion the transition to the ms = level is active. More de- 
tailed analysis of this mechanism is given in IV A 31 where 
we associate this with a magnetic dipole transition. 

Within the error limits the g- factors of the triplet state 
are the same for the two in-plane field orientations g a — 
1.97 ± 0.02, g b = 1.975 ± 0.004. The fit of the resonance 
frequencies of the singlet to triplet transition for the third 
field orientation, Bo || c, gives g c = 1.90 ± 0.03. 

In zero field the triplet levels ms — —1,0, and 1 are 
degenerate; it is best seen when comparing the zero field 
line positions in Fig lllf c) and (d). Determining the size 
of the zero field splitting is limited by the linewidth. We 
can say that the zero field splitting of the triplet levels is 
less than half of the linewidth, FWHM/2=0.25cm -1 . 
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FIG. 14: The magnetic field dependence of the singlet to 
triplet transition line area at 4.4 K in two polarizations, Ei || c 
(empty symbols) and Ei || b (filled triangles). The line areas 
are plotted in the negative field direction for the transition to 
the ms = — 1 triplet level and for the transition to ms ~ 1 
in the positive field direction. The zero field data points on 
the graph are one half of the measured area, (a): Bo || a and 
Bo || b. (b): Bo || c. The dashed line in (a) is a guide for 
the eye. The dashed line in (b), units on the right axis, is 
the 68 cm -1 absorption line shape with the background sub- 
tracted and the energy units converted into magnetic field 
units using the triplet state g- factor, g c = 1.90. The error 
bars shown only in (b) apply to data points in both panels. 
A solid line in (b) is a fit to the dynamic DM mechanism in- 
duced by the 68 cm -1 phonon. The oscillator strength of the 
phonon and the singlet to triplet transition in zero field are 
400 cm -2 and 0.13 cm -2 , respectively. The only fit parameter 
is the dynamic DM interaction, gDq = 0.13 cm -1 . 



V. DISCUSSION 

A. Triplet state and Dzyaloshinskii-Moriya 
interaction 

In this section we analyze three different contributions 
to the singlet to triplet optical absorption observed exper- 
imentally in a'-NaX^Os. Two of them fall into the same 
category, the electric dipole active dynamic DM mech- 



anism. The third contribution is probably due to the 
static DM interaction mechanism where magnetic dipole 
transitions are active, but cannot be explained with a 
single dimer model. 

We rule out a possible singlet to triplet absorption 
mechanism based on a staggered g-factor. The staggered 
g-factor mechanism requires that the principal axes of 
the g-tensors of a pair of spins with anisotropic g-factors 
must not coincide. Let us consider for an example the 
spin Si with its principal axes rotated from the crys- 
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tal a-axis by an angle 9 and the spin S2 by an angle 
— 9 and B || a. Let the hamiltonian contain isotropic 
exchange interaction JSi • S2 and Zeeman interaction 
Bo • (gi • Si + §2 • S2). Matrix elements between the 
singlet and triplet states equal to ±B (g a — 9b) sin 9 cos 9 
appear. It is important that the singlet-triplet mixing is 
proportional to Bq. First, there is no mixing in zero field 
(So = 0) and the staggered g-factor mechanism is turned 
off. Second, the transition probabilities increase as (B ) 2 
if the magnetic held is small compared to the separation 
of the singlet and triplet states. In the experiment we 
observe a singlet to triplet transition in zero field, and 
also the observed magnetic field dependence is different 
from that of the staggered g-factor mechanism. There- 
fore this mechanism does not apply to the singlet-triplet 
transitions in c/-NaV205. 



1. Dynamic Dzyaloshinskii - Moriya: Ei ||a 

The strongest singlet to triplet absorption is observed 
in Ei || a polarization. There is no resonant enhance- 
ment in the magnetic field dependence of the line area 
as seen in Fig.^] and any of the optically active excita- 
tions in Ei || a polarization is a candidate that can create 
the dynamic DM interaction. Nevertheless we can make 
some choices. The 91 cm -1 a-axis phonon is not active 
since no enhancement is observed when the upper branch 
of the triplet resonance crosses the phonon frequency at 
28 T fFig. 1121) although there is an interaction between 
this phonon and the continuum of magnetic excitations 
as is manifested by the Fano line shape of the phonon 
line above T c (Fig.[5J. 

In Fig. E| we have plotted two fit curves based on 
perturbation calculation results (Eq.-s I14|l . In one 
case the phonon frequency was fixed to 518 cm -1 that is 
the strongest a-axis optical phonon^i. In the other case 
the phonon frequency was a fitting parameter giving us 
(180± 10) cm" 1 . This value represents the lowest bound- 
ary for the frequency of the DM phonon. Phonons with 
lower frequencies would give a too steep magnetic field 
dependence of the singlet to triplet transition probabil- 
ity. It is likely that the lowest boundary of the phonon 
frequency has been underestimated since the three high 
field data points above 30 T influence the fit by lowering 
the phonon frequency. Additional measurements above 
33 T are required to clarify this intensity enhancement. 
Based on these fits the strength of the dynamic DM in- 
teraction qDQ can be obtained. The closest low T phase 
infrared active a-axis phonon in frequency to 180 cm -1 
is the 199 cm -1 zone folded phonon^I. The plasma fre- 
quencies fi p of the 199 cm -1 and 518 cm -1 phonons are 
48 cm -1 (Refi2£) and 853 cm -1 (RefiSi), respectively. We 
convert the plasma frequency into the line area in ab- 
sorbance units, J a(u>)dw, using f a(u>)du) — ir 2 flp/n a . 
From the fit results we calculate qDQ = 5 cm -1 for the 
199 cm -1 phonon and qDQ = 0.9 cm -1 for the 518 cm -1 
phonon. 
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FIG. 15: Cartoons of superexchange paths in the ab-plane 
(panels (a) and (b)) with one ladder shown and in the ac- 
plane ((c) and (d)) with two ladders shown. For illustrative 
purposes it is assumed that the spin is located on the rung 
oxygen. Oxygens are shown by filled circles and vanadium 
atoms by open circles. Block arrows show the displacement of 
atoms due to a phonon. (a) superexchange paths between two 
spins over on-leg oxygens within the same ladder in the zig- 
zag ordered phase; (b) same superexchange paths when an a- 
axis phonon is involved; (c) superexchange paths between two 
spins in the neighboring planes over apical oxygens displaced 
by a c-axis phonon when spins are in-line along c-axis or (d) 
zig-zag. 

The orientation of the DM vector is determined by the 
vector product of the vectors connecting two spins over 
the superexchange pat h 34 i 35 

Pi = [R-SlO, x R 0l ,s 2 ] (16) 

In the zigzag ordered low T phase two superexchange 
paths exist, one over on-leg oxygen Oi and the other 
over O2, shown in Fig |15f a). The resultant DM vector 
is zero since the two DM vectors, pointing in c-direction, 
cancel each other. The zero length of the DM vector 
follows from the general arguments of symmetry too as in 
this particular case there exists a local center of inversion 
located between two V-O-V rungs. The 518 cm -1 a-axis 
phonon displaces on-rung oxygens along the rung (on- 
rung V-O-V stretching mode)2ii, Fi sHSf bh Because of 
the phonon the two superexchange paths are not "equal" 
any more and the resulting DM vector points along the c- 
axis. The orientation of the dynamic DM vector, Dq || c, 
is consistent with the selection rules for the dynamic DM 
interaction observed experimentally. The transition to 
the |t 0) is observed when B || Dq, Fig.HlTdh The 
transitions to the |t_0) and |t+0) states are observed 
when B J_ D Q , Fig.EHa), (b), (c). 

There is another phonon with an oscillator strength 
larger than that of the 518 cm -1 a-axis phonon. That 
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is the 582 cm -1 b-axis phonon that stretches on- leg V- 
O-V bonds 31-33 . In the high T phase where electron 
charge (and spin) is rung-centered, such distortion does 
not produce any DM interaction from the principles of 
symmetry. In the low T phase the charges are ordered 
in a zigzag pattern. The on-leg V-O-V bond stretching 
phonon lowers the symmetry and creates the DM inter- 
action along the c-axis. We do not observe a singlet to 
triplet absorption with an electric dipole moment along 
the b-axis in the experiment. There are two possibili- 
ties why the Ei || b absorption is not observed. First, 
there is no zigzag charge order. Second, there is a zigzag 
charge order and although the dynamic DM is allowed 
by symmetry, the actual value of qDq is small and the 
absorption cannot be observed in our experiment. Since 
there is a mounting evidence in the favor of a low T zigzag 
charge ordenZ&il we consider the second case likely. 



2. Dynamic Dzyaloshinskii - Moriya: Ei ||c 

The singlet to triplet absorption with the electric 
dipole moment along the c-axis is due to the dynamic 
DM interaction. In this particular case the dynamic DM 
mechanism is brought in by the 68 cm" 1 optically active 
c-axis phonon. The enhancement of the singlet-triplet 
absorption close to the 68 cm" 1 line is present if the mag- 
netic field is either parallel to a- or b-axis, but missing if 
the field is along the c-axis. In Section lY Cl we present fur- 
ther arguments supporting the assignment of the 68 cm -1 
resonance to a phonon and not to a magnetic excitation. 

In Fig.^l the fit of the singlet to triplet transition 
intensities to the dynamic DM absorption mechanism 
is shown. The input parameters are the resonance fre- 
quency of the phonon, uj p = 68 cm -1 , and the frequency 
of the singlet to triplet transition as the function of mag- 
netic field, lot ± = AzLg/iBBo, where A = 65.4cm" 1 . We 
estimated the low temperature oscillator strength of the 
phonon from the T dependence presented in Fig. llUI and 
got Hp = 400 cm" 2 . The ratio of the singlet to triplet ab- 
sorption oscillator strength f2| T to the oscillator strength 
of the phonon, ft 2 ,, is ^| T /ft 2 ~ 0.13/400 = 3.3 x 10" 4 . 
The only fit parameter is the strength of the dynamic 
DM interaction, qDQ = 0.13 cm" 1 . 

We used Dq || b in our fit because the splitting of 
the triplet in the magnetic field is observed if Bo || a or 
Bo || c fFig. H4fl . According to the selection rules for the 
dynamic DM only the transitions to the triplet states 
with ms — 1 and m§ = — 1 are observed when the mag- 
netic field is perpendicular to the dynamic DM vector 
Dq _L Bo. The selection rules for the dynamic DM if 
Bo || Dq I) b allow only transitions to the mj = level 
that does not shift with the magnetic field and we do not 
see it in the differential absorption spectra that are taken 
in different magnetic fields. The background intensities 
that do not depend on the magnetic field in figure H4T a1 
are analyzed in the next section (|V A 31) . 

Which lattice deformations along the c-axis (electric 



dipole!) will give the dynamic DM interaction in the b- 
axis direction? It turns out that we have to consider 
inter-plane interactions. For the beginning let us con- 
sider two spins on the neighboring rungs, as shown in 
Fig |15f ah where optical c-axis phonons have the out of 
(ab) plane anti phase movements of oxygen and vana- 
dium atoms. In the low T zigzag ordered phase the dy- 
namic DM vector will have components along both, a- 
and b-axis. The b-axis component of the dynamic DM 
turns to zero for a vanishing zigzag order, i.e. in the limit 
of rung-centered spin distribution. The a-axis component 
stays nonzero. To get a dynamic DM exclusively along 
the b-axis we have to consider interactions between the 
spins in the neighboring planes. Here the path for the 
DM interaction between spins in the neighboring planes 
goes over the apical oxygens. Two arrangements along 
the c-axis are possible in the zigzag ordered phase, shown 
in Fig ll5f c).fdh in-line or zigzag. In both arrangements 
the displacement of the apical oxygens in the c-direction 
will create a dynamic DM along the b-axis. As one can 
see not only the dynamic, but also the static DM in the 
b-direction is allowed by the symmetry. 



In the experiment we do not see magnetic dipole active 
optical transitions caused by the static DM interaction. 
To compare the strength of magnetic and electric dipole 
transitions we must know the magnitude of the static DM 
interaction. The magnitude of the DM interaction was 
estimated by MoriyaSi, D ss (\g — g e \/g e )J, where J is 
isotropic exchange interaction and g e — 2 is the free elec- 
tron g-factor. We use g = 1.90 and J = 60meV (from 
Refi 1 *) and get D = 25 cm" 1 . The transition intensity 
I is proportional to (|F| /(!>£) 2 where |V| is the matrix 
element of the interaction (D or <7-Dq) between the two 
states and S£ is their energy separation. For the static 
DM effect S£ is the singlet-triplet gap, S£ — 65.4 cm" 1 ; 
for the dynamic DM it is the energy difference between 
the phonon energy and the triplet level, 8£ — 68 — 65.4 = 
2.6 cm -1 . The intensity of the magnetic dipole transi- 
tion is weaker than the electric dipole transition by the 
factor of exj 1 = 137 2 , where a/ is the fine structure 
constant (see Refi2£ p. 171). If we use D = 25cm" 1 and 
qDQ = 0.13 cm -1 we get that the intensity due to the dy- 
namic mechanism is (0.13 x 65.4) 2 /(a/25 x 2.6) 2 « 300 
times larger than the intensity due to the static mecha- 
nism. Optical transitions due to the static DM interac- 
tion are suppressed with respect to the transitions caused 
by the dynamic DM because electric dipole transitions 
are stronger than magnetic dipole transitions. 



Our conclusion is that the dynamic DM interaction 
along the b-axis, qDQ = 0.13 cm" 1 , is between the spins 
on the ladders of the neighboring planes. The 68 cm -1 
c-axis phonon that creates the dynamic DM interaction 
involves the displacement of apical oxygens. 
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3. Third optical singlet-triplet absorption mechanism 

There is a third mechanism for the optical triplet ab- 
sorption that is responsible for the magnetic field inde- 
pendent intensities of the transitions to the ms = -l and 
?7i5 = 1 triplet levels, shown in Fig. ll4f a) by solid trian- 
gles and empty circles. The same mechanism contributes 
together with the enhanced part discussed in the previ- 
ous section to the intensity plotted with empty triangles 
on the same graph. 

We could assume that the third mechanism is also an 
electric dipole absorption mechanism as the two mecha- 
nisms ascribed to the singlet-triplet absorption in Ei || a 
and Ei || c polarizations. Since this absorption is present 
in both polarizations, Ei || b and Ei || c one has to as- 
sume that the electric dipole moment is either in the (bc)- 
plane or, just by coincidence, two electric dipole mecha- 
nisms, one polarized along the b-axis and the other polar- 
ized along the c-axis, give the same intensities. The case 
that the optical phonon responsible for the dynamic DM 
effect has a dipole moment in the (bc)-plane contradicts 
with the data available on phonons and with the crystal 
symmetry. The second case of coinciding intensities is 
ruled out by the selection rules if applied to the full data 
set presented in Fig. ll4f a) and (b). 

The third mechanism could be a magnetic dipole DM 
mechanism. By applying the selection rules to the data 
we should be able to determine the orientation of the DM 
vector. According to the theory fsection fll Aj l transitions 
to the 777,5 — — 1 an( l m s — 1 states have constant and 
equal intensities when Bo || D. This condition is satis- 
fied by two data sets, solid triangles (Hi || c) and empty 
triangles (Hi || b) in Fig. ll4f a). measured with Bo || a. 
We have D || a, where if = if and x = b and y = c 
(see Fig.^l. However the set represented by the circles 
in Fig. ll4f a) is the I~ intensity in the perpendicular con- 
figuration, Bo J- D (in our notation D || y || a), which 
should have zero intensity according to the theory (see 
Fig. |2| • Also, in the Bo || Hi || c configuration if should 
have non-zero intensities. In the experiment, filled trian- 
gles in Fig. ll4f b). no intensity is observed contrary to the 
theory. Therefore the assumption D |] a is not consistent 
with the full data set. Also, we can prove that neither 
D || b nor D || c is fully consistent with the experiment. 

Our conclusion is that the third optical triplet absorp- 
tion mechanism is a magnetic dipole transition, but can- 
not be explained by an isolated dimer model with DM 
interactions. Apparently, a more elaborate model includ- 
ing the inter-dimer DM interactions, is necessary. 

B. Spin gap and phase transition 

We studied the effect of temperature and magnetic field 
on the singlet to triplet transition. The results arc shown 
with circles in Fig.[S] We observe that the intensity of the 
singlet to triplet transition (Fig.|SJa)) follows the inten- 
sity of the X-ray diffraction peak2& reported in the same 
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FIG. 16: Relative deviation (E-E fit )/E of S to T_i (plotted 
at negative field values) and S to T+i (plotted at positive field 
values) transition energy E from the linear fit -B/it- 

figure. Gaulin et al. have shown that the intensity of the 
X-ray diffraction peak can be fitted with a single func- 
tion I = t° over the reduced temperature, t = 1 — T/T c , 
in the range from 6 x 1CP 3 to 2 x 10 -2 with the critical 
exponent (3 = 0.18. Our data taken below 30 K is above 
1.1 x 10 _1 in the reduced temperature scale, which is 
unsuitable for the determination of the critical exponent. 

The singlet-triplet splitting (spin gap), shown in 
Fig-E^b) , is more rigid than the singlet to triplet tran- 
sition probability or the X-ray scattering intensity. The 
rigidity of the spin gap has been confirmed earlier by 
INS measurements^ and by high field electron spin reso- 
nance measurements, although not in the zero magnetic 
fielcU^. The far-infrared S to T gap is in agreement 
with the INS gap also shown in Fig.[HJb). In the ul- 
trasonic experiment^ the measurements were extended 
close enough to T c that one is able to determine the crit- 
ical exponent for the spin gap. The temperature depen- 
dence of the spin gap, measured indirectly by the ul- 
trasonic probe, gives the critical exponent (3 — 0.34 be- 
low the reduced temperature 10~ 2 . Although far- infrared 
measurements of the gap are direct, reasonable data can 
be obtained only too far from T c as the absorption line 
broadens and gets weak and therefore the fit, inset to 
Fig.[S{b), gives us an exceptionally low critical exponent, 
13 = 0.039 ±0.002. 

The magnetic field dependence of the spin gap and the 
triplet state g-factor have been measured up to 33 T for 
one field orientation, Bo || b, at 4.7 K. The line posi- 
tions of the transitions from the singlet ground state to 
the triplet states with ms = ±l were fitted with a linear 
function Efu = A± g/isBo- The reduced residual of the 
linear fit is less than 0.0012 of the transition energy, as 
shown in Fig.^J This means that the triplet state g- 
factor is not renormalized by the magnetic field as high 
as 33 T and the spin gap is independent of the field up to 
33 T. The independence of the gap value A on the field at 
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low temperature is not surprising 40 . The role of the mag- 
netic field is to couple to thermally excited quasiparticles. 
At low temperature the number of quasiparticles is small 
and hence the effect of the magnetic field is negligible. 

In conclusion, the temperature dependence of the in- 
tensity of the singlet to triplet transition agrees with 
the X-ray scattering intensity temperature dependence. 
Also, the singlet-triplet splitting at 8.13 meV has the 
same temperature dependence as the singlet-triplet split- 
ting of the second excitation branch at 9.8 meV measured 
by the INS. The 8.13 meV spin gap is not altered by the 
magnetic field at least up to 33 T. 



C. Phonons and phase transition 



1. c-axis phonons 



In Section lY A 2l while calculating the dynamic DM in- 
teraction for Ei || c we assumed that the electric dipole 
moment of the c-axis polarized optical singlet to triplet 
transition comes from the interaction between the spin 
system and the 68 cm -1 optical c-axis phonon. An al- 
ternative would be that the 68 cm -1 resonance is not a 
phonon but a singlet electronic excitation. Here we ana- 
lyze existing data and show that the data is not in contra- 
diction with the assumption that the 68 cm -1 resonance 
is a phonon mode. 

There is a series of infrared-active c-axis modes ap- 
pearing below 34K as shown in Fig.[S]and found by other 
groups 32-41 . The two low frequency modes, at 68 and 
106 cm" 1 , are the strongest. The temperature depen- 
dence of their resonance frequencies plotted in Fig. lT7T a') 
is an order of magnitude larger as compared to the a- 
and b-polarized modes plotted in Fig. ll7f b1. Modes at 
frequencies close to 68 and 106 cm" 1 have been found by 
Raman spectroscop y 4 ^ 43 . Sample dependence of mode 
frequencies by few wavenumbers has been observed both 
in infrared 3 ^ and RamanS. It is known that Na defi- 
ciency of a'-Na^^^Os affects T c and in samples where 
x = 0.97 the phase transition is suppressed completely 44 . 
A Na doping dependence study^ has shown that the fre- 
quency variation of Raman modes from sample to sam- 
ple is related to the Na content of the sample. Raman 
measurements^ on the sample from the same batch as 
the one studied here gave values for the mode frequen- 
cies 65.9 and 105 cm" 1 at 4.4 K. This means that Raman 
and infrared spectroscopy observe different modes since 
the frequencies of infrared modes are 68 and 106 cm" 1 . 
Neither of the Raman modes^ 2 ^ 3 - nor the infrared modes 
(this work and Ref4i) split in the magnetic field. The 
comparison of the T dependence of the normalized mode 
frequency has shown that the T dependence of the 
65.9 cm" 1 Raman mode follows the T dependence of the 
68 cm -1 infrared mode and the 105 cm" 1 Raman mode 
follows the T dependence of the 106 cm" 1 infrared mode. 
Therefore, although the infrared and the Raman modes 
do not have the same frequencies, it is likely that the 
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FIG. 17: Temperature dependence of normalized phonon fre- 
quencies, (a) - b-axis multiplet and c-axis phonons, (b) - a- 
and b-axis phonons. 



origin of the infrared and Raman modes near 68 and 
106 cm -1 is the same, lattice modes or electronic excita- 
tions. 

The origin of the low frequency infrared and Raman 
modes has been under the debate. It has been concluded 
that the Raman modes are not lattice vibrationals^ 3 .. 
In RefiS the conclusion was that the infrared modes at 
68 and 106 cm -1 are zone-folded c-axis phonons and it 
was speculated that their peculiar T dependence com- 
pared to other zone-folded phonon modes is caused by 
the interaction of phonons with charge and spin degrees 
of freedom. Indeed, the 68 cm" 1 mode and the singlet to 
triplet excitation at 65.4 cm -1 have similar T dependen- 
cies of their infrared absorption line parameters as shown 
in Fig.|S| Assuming that 68 and 106 cm -1 excitations are 
zone folded phonons we would have to explain why their 
frequencies depend on Na deficiency and why their fre- 
quencies have a different temperature dependence than 
the other zone-folded modes have. We argue that both 
points can be explained if charge correlations develop in- 
side ladder planes prior to the transition to the low T 
phase and at the phase transition three-dimensional cor- 
relations build up between the planes. 

Charge ordering within ladder planes precedes the lat- 
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tice distortion and the opening of the spin gap as ev- 
idenced by the vanadium Knight shift and the sodium 
quadrupolar and Knight shifts^. The presence of two- 
dimensional charge correlations in the ladder planes 
above 35 K is also supported by the X-ray diffraction 
measurements^. In a'-NaX^Os there are modulated lad- 
ders with a zigzag charge order within one plan o^i?! 1 ? . 
It is possible to construct four different planes with the 
zigzag charge order and when stacked in certain sequence 
along the c-axis a unit cell is formed compatible with the 
observed X-ray structur e^ 11 . By applying pressure the 
critical temperature is reduced and it has been found that 
the soft axis is the c-axis^^ 6 .. A "devil's staircase" like 
sequence of phase transitions to the phases with unit cells 
incorporating more than four planes in the c-direction 
takes place under pressure larger than 0.5 GPa£L. There- 
fore the interactions between layers are important in the 
formation of three dimensional correlations. 

Interactions between ladder planes are important for 
those modes, which are zone-folded along the c-axis. The 
primary candidates are the modes with anti-phase move- 
ments of atoms along the c-axis in the neighboring planes, 
which are naturally the modes that are polarized along 
the c-axis. The c-axis modes have the largest frequency 
shift close to the phase transition point. The exception 
is the b-polarized multiplet (Fig.^J . It is the only non-c- 
axis low frequency mode, which has a relative change of 
frequency with temperature (Fig. ll7f all that is compa- 
rable to the c-axis modes. Considering the low frequency 
of the multiplet it is likely that the seven lines are a re- 
sult of the folding of a b-polarized acoustical, not optical 
phonon branch. To get seven optical phonon branches 
out of one phonon branch a c-axis folding is required in 
addition to a- and b-axis foldings. Therefore this multi- 
plet involves relative movements of atoms in neighboring 
layers and is a subject to interlayer couplings. 

Our conclusion is that the c-axis phonon modes are 
most susceptible to three dimensional correlations be- 
tween ladder planes and therefore their T dependence 
is different from other low frequency phonon modes. We 
assign the 68 cm- 1 (and 106 cm -1 ) infrared active mode 
and the 66 cm -1 (and 105 cm -1 ) Raman active modes to 
the zone- folded phonons. 



2. Polarization dependence of line intensities 

In NaV2C>5 the a- and b-axes are not equivalent and 
one would expect the selection rules to apply so that 
the phonon active in the Ei || a absorption is missing in 
Ei || b spectrum and vice versa. In other words, phonons 
polarized along a-axis should have different frequencies 
than phonons polarized along the b-axis. 

This rule applies well at 40 K (Fig. The strong 
180 cm -1 phonon is missing in the Ei || a spectrum and 
the structure at 140 cm -1 is missing in the Ei || b spec- 
trum. The 90.7 cm -1 line is present only in the Ei || a 
spectrum. 
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FIG. 18: Polarization dependence of low frequency phonon 
absorption line intensities, (a) - absorption spectra in Ei ||a 
(shifted upwards by 150 cm ) and Ei ||b polarizations, (b) - 
absorption line areas as a function of polarization angle in the 
(ab)-plane; r\ = corresponds to Ei || b. Inset to (a) shows 
the orientation of electric dipoles rotated from the b-axis by 
an angle <p and — tp. 



At 4.4 K this selection rule does not apply to all absorp- 
tion lines. The 101.4 and 111.7, 126.7 and 127.5 cm- 1 
phonons (Fig.-s ITHT ah 151 have absorption in both polar- 
izations. The leakage of the wrong polarization through 
the polarizer is ruled out. For example, the 101.7cm -1 
line present in Ei || a spectrum is missing in Ei || b spec- 
trum (Fig. 118b) and the multiplet at 30 cm -1 is missing 
in Ei || a spectrum (Fig. HI . 

We measured the polarization dependence of line in- 



TABLE II: Rotation angle ip calculated from the two dipole 
model, (fa (tfb) is the rotation angle of the phonon dipole 
moment from the crystal a-axis (b-axis). 
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tensities in the (ab)-plane. The areas of absorption lines 
between 100 and 130 cm -1 are plotted in Fig. I18f b) as 
a function of the angle 7/ between the polarizer and 
the b-axis. The lines at 91 cm -1 were not included 
since they are two different lines, one at 91.2 cm -1 in 
the a-polarization and the other at 91.3 cm -1 in the b- 
polarization. As the 101.7cm -1 line has zero intensity at 
r] = 0°, Ei || b it belongs to the phonon modes represent- 
ing the averaged crystal symmetry. It is important that 
the intensity of the 101.4cm -1 , 111.7cm -1 , 126.7cm -1 
and 127.5 cm -1 lines never falls to zero between r\ = 0° 
and 90°. A similar effect was found by Damascelli et 
al£L for several infrared active phonons at higher fre- 
quencies. The simplest explanation would be that the 
zone-folded a- and b-polarizcd phonons are pair wise de- 
generate, but then it is hard to explain why there is a 
constant intensity ratio of at least four pairs of low fre- 
quency and of one high frequency (718 cm -1 ) pair of a- 
and b-axis phonons, as is shown below. 

We assume, as was done in Refill, that there are two 
regions in the crystal with a symmetry different from the 
averaged crystal symmetry. In one region the electric 
dipole moment of a phonon, d, is rotated from the crys- 
tal b-axis by an angle tp and in the other region by — tp as 
shown in the inset to Fig ll8f a). The optical conductivity 
in the first region is (Ji{<p) = aE 2 d 2 cos 2 (tp — rj), where 
a is a constant independent of angular parameters. 77 is 
the angle between the b-axis and the electric field vec- 
tor of light Ei . The conductivity in the second region is 
(Ti (— ip) = aE\d 2 cos 2 (tp+rj). Since the phonons are inde- 
pendent in the two regions the total observed conductiv- 
ity is (Ti = <Ti(<p>) + ci(— tp) — lb cos 2 r\ + I a sin 2 77, where 
lb = 2aE\d 2 cos 2 tp and I a — 2aE 2 d 2 sin 2 tp are the oscil- 
lator strengths observed in the experiment in the b- and 
a-polarizations respectively. The solid lines in Fig. l^f b - ) 
are the fits of the integrated absorption line areas where 
Eg . 1151 has been used to convert a(uj) into ui(uj). The 
results are given in Table [HI Note that two of the ab- 
sorption lines belong to the dipoles rotated from the a- 
axis and two to the dipoles rotated by approximately the 
same angle from the b-axis. The averaged value of the 
rotation of the dipole moment is (p = 34° ±3°. A simi- 
lar result, ip = 39°, we obtain for the 718cm -1 phonon 
using the data from the Refi^i where the given oscillator 
strength in the Ei || a polarization is 0.021 and in the 
Ei || b polarization is 0.014. 

The two regions with a symmetry that is different from 
the crystal symmetry can be associated with two types 
of ladder planes with a zigzag charge order. In Ref^i it 
was assumed that the crystal is split into domains with 
two different diagonal charge patterns. Considering the 
recent X-ray studie a 8 ' 11 it is more likely that there are 
no domains and two types of planes exist with different 
diagonal charge pattern instead. The two different charge 
configurations come from the way the two neighboring 
zigzag charge ordered ladders are positioned with respect 
to each other within the plane. A shift of every second 
ladder by one half of the superlattice constant in the b- 



direction creates two different structures where charged 
stripes of v +4 ' 5-<5c / 2 run diagonal from left to the right or 
from right to the left in the (ab)-planeii. This diagonal 
charge order within a single ladder plane determines the 
orientation of the phonon dipole moments. It is possible 
to estimate the charge offset x = 1(1 + S c )/2 from the 
center of the rung knowing the tilt angle of the dipole 
moments tp. Here I is the length of the rung and the 
charge transfer factor S c is associated with the formal 
valence of vanadium ions in the low T zigzag ordered 
phase, V+ 4 5- ^/ 2 and V +4 - 5+<5 <=/ 2 . The charge transfer 
factor depends only on the angle tp, S c ~ tan tp, since 
the rung length is approximately equal to the distance 
between the nearest neighbor V atoms along the leg. We 
get S c = 0.67 ± 0.07 using tp = Cp. From the analysis of 
INS data the authors of Refpi^ come to a similar value, 
S c = 0.6. 

Our study of phonon modes in the low temperature 
phase supports the view that the symmetry of the indi- 
vidual ladder planes is lower than the averaged crystal 
symmetry. In addition it has been found by sound veloc- 
ity measurements^ that the cee shear mode couples to 
the pre-transitional charge fluctuations of B\ g symmetry, 
which correspond to the static zigzag charge order in the 
low T phase. There are phonon modes in a'-NaV205 
in the low T phase where the normal coordinates are 
confined into single ladder planes. These modes show 
the symmetry of an individual plane determined by the 
zigzag charge ordering. For other modes the movement 
of atoms is correlated between the neighboring planes or 
they are insensitive to the zigzag charge order and there- 
fore they reflect the averaged crystal symmetry. 

From the analysis of the infrared spectra we conclude 
that two types of ladder planes exist with the zigzag 
charge order where the charged stripes are aligned ap- 
proximately in [110] or in [110] directions. The formal 
charge of the vanadium atoms forming the zigzag pat- 
tern is +4.17 ± 0.04 and +4.83 ± 0.04. 



D. Continuum of excitations and Fano resonances 

A broad absorption continuum is observed in Ei || a po- 
larization above T c that starts gradually below 20 cm -1 
(Fig. Ela)) and extends up to 400cm" 1 (RefA). The 
change of the continuum with temperature is correlated 
with the phase transition at 34 K as demonstrated in 
Fig. 1191 where normalized optical conductivity is plotted 
for three different frequencies. The continuum absorp- 
tion at 4.4 K shows a threshold at 130cm -1 in Ei || a 
polarization while at lower frequencies it is similar to the 
absorption in Ei || b polarization fFig. I20fl . It is natural 
to associate the absorption continuum with an optical 
excitation of two triplets since the lowest in energy one 
triplet excitation is at 65.4cm -1 . 

There are two derivative- like phonon absorption lines, 
one at 91.2 and the other at 140 cm -1 , in the Ei || a ab- 
sorption spectrum above T c in Fig.|SJa). This line shape, 
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FIG. 19: Temperature dependence of the normalized opti- 
cal conductivity in Ei || a polarization uin = [o\ a (T) — 
aia(4.4A")][cria(40A')-cria(4.4A")]- 1 at 60 (open circles), 108 
(crosses), and 160cm -1 (filled circles). 
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FIG. 20: Real parts of optical conductivities <j\ a (shaded and 
thin line) in Ei || a polarization and au (thick line, Ei | b) 
at 4.4 K. Arrows point to the singlet to triplet resonance at 
65.4 cm -1 and to the onset of the absorption continuum at 
132 cm -1 in Ei || a polarization. 



known as the Fano resonance, comes from the interaction 
of a discrete level and a continuum of state d 49 ! 50 . The 
Fano lineshape^ - does not describe a normal absorption- 
like lineshape. To account for an arbitrary line shape, in- 
cluding the lorcntzian line, we use an empirical formula^ 
for the dielectric function e(w) where the asymmetry is 
described by the phase 9: 



n 2 p exp(-i0) 

LU 2 — LU 2 — ILJ^j 



(17) 



dp is the plasma frequency, e m is the background di- 
electric constant, and u p is the resonance frequency of 
the phonon, and 7 is the full width at half maximum 
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FIG. 21: Temperature dependence of the 91cm~ phonon 
absorption line in Ei || a polarization. Spectra have been 
off-set by 50 cm -1 in vertical direction. 



(FWHM) of the lorentzian line at 9 = 0. 

We did a T dependence study of the 91 cm -1 line and 
the evolution of this line with T is shown in Fig. 1211 
The phonon line at 140 cm -1 is on a steeply rising back- 
ground. Since the frequency dependence of the absorp- 
tion background is not known it is difficult to subtract 
the background reliably and we did not attempt to eval- 
uate the T dependence of the 140 cm -1 line parameters. 
The absorption lines were fit with Eq.^l The relation 
between the real part of the conductivity and dielectric 
function is (Ji(lu) — to Im[e(aj)] /60, where <j\ is in units 
il~ 1 cm _1 and Eg .11 51 was used to calculate o\ from the 
measured absorption spectra. The fit parameters are 
plotted in Fig. [23 The phase 9 (Fig. ICTdD changes 
from 3ir/2 at 35 K to n at 33 K where the line has the 
shape of an antiresonance (Fig. l21f) . At 31 K the line has 
a derivative-like lineshape again, but with a phase 7r/2. 
Also the phonon resonance frequency uj p stays fairly con- 
stant at higher temperatures and changes abruptly at 
35 K. The line area S and 7 have a smooth T depen- 
dence around T c — 35K. At 4.4 K 9 = and the normal 
lineshape is recovered. The line width 7 (Fig. 122( a)) goes 
down with T and is limited by the 0.4 cm -1 instrumental 
resolution for this T dependence study. We know from 
the higher resolution measurement that the line width is 
0.2 cm -1 or less at 4.4 K. The line area (Fig.^c)) has a 
different T dependence than other parameters, having a 
minimum at approximately 32 K. S changes substantially 
even above T c , by a factor of 2 from 35 to 40 K. 

In this paper we will not present a theory covering the 
optical conductivity of a system where phonons interact 
with a two-particle continuum of magnetic excitations. 
A proper theory must account for a microscopic mecha- 
nism responsible for the optical absorption continuum 5 J 
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FIG. 22: Temperature dependence of the Fano fit parameters 
of the 91 cm -1 phonon line: (a) - full width at half maximum 
7, (b) - phonon resonance frequency lj p , (c) - integrated line 
area 5, (d) - asymmetry parameter 6. Open circles - 120 /im 
thick sample; filled circles - 40 /jm thick sample. 



Nevertheless, some observations can be made based on 
the empirical fit of the phonon line shape. The phase 9 
that is related to the Fano parameter q, g" 1 oc tan(0/2), 
depends on the strength of the interaction between the 
phonon and the magnetic system. The T dependence of 
shows that the spin-phonon interaction weakens as T 
is lowered below 35 K. The other reason why the normal 
lineshape of the 91 cm -1 phonon at low T is recovered 
could be vanishing of the two particle absorption con- 
tinuum below 130 cm -1 . The direct evidence that the 
spin-phonon interaction is switched off or is very weak at 
low temperatures is the lorentzian shape of the 140 cm -1 
phonon line. 

In conclusion, our data shows that there is an absorp- 
tion continuum in the a-axis polarized optical absorp- 
tion that develops a gap at low T. This gap, 130cm" 1 , 
is equal to twice the singlet-triplet excitation energy, 
65.4 cm -1 . Therefore the absorption continuum can be 
assigned to an absorption of a photon with a simulta- 
neous creation of two magnetic (spin) excitations. Two 



low frequency a-axis polarized optical phonons, at 91 and 
140 cm -1 interact with the magnetic system as evidenced 
by their derivative-like line shape in the high T phase. 



VI. CONCLUSIONS 

Using far-infrared spectroscopy we have probed spin 
excitations and phonons in the quarter-filled spin lad- 
der compound a'-NaV205. The interaction between the 
spins and the phonons is observed in the gapped spin 
state below 34 K and above, in the paramagnetic phase. 

The zigzag charge order within the ladders in the 
gapped state is in accordance with the polarization de- 
pendence of several infrared active zone-folded optical 
phonons. In the high T phase two a-axis optical phonons 
interact with the continuum of excitations as manifested 
by the Fano line-shape of phonon lines. At low T the 
continuum absorption is gapped with a threshold energy 
130cm" 1 , twice the singlet-triplet gap, and the interac- 
tion between the phonons and the continuum of excita- 
tions is turned off. 

The strength of the singlet to triplet absorption at 
65.4 cm" 1 is strongly anisotropic. Absorption is strongest 
when the electric field of the incident light is polarized 
along the ladder rungs (Ei |j a). In this polarization 
the strength of the singlet to triplet absorption has a 
weak magnetic field dependence up to 28 T. In Ei || c 
polarization a strong magnetic field dependence of the 
singlet-triplet absorption is observed. This field depen- 
dence is due to the dynamic DM interaction created by 
the 68 cm -1 c-axis optical phonon. In the case of the dy- 
namic DM absorption mechanism the singlet-triplet ab- 
sorption is electric dipole active and the polarization of 
the transition is determined by the polarization of the op- 
tical phonon creating the dynamic DM interaction by the 
lattice deformation. Using the presented theory we cal- 
culated the strength of the dynamic DM interaction. The 
dynamic DM interaction, qDQ = 0.13cm" 1 , created by 
the 68 cm -1 c-axis optical phonon is inter-ladder between 
the spins in the neighboring planes and points along the 
b-axis. We assign the strong a-axis polarized absorption 
also to a dynamic DM effect, in this case due to the lattice 
deformation caused by one of the high frequency a-axis 
optical phonons. This dynamic DM interaction is intra- 
laddcr and is along the c-axis. Above 28 T an increase in 
the S to T + a-axis absorption intensity is observed. The 
origin of the mechanism responsible for that increase is 
not clear and additional measurements above 33 T are 
required. 

To summarize, the optical singlet to triplet transition 
in a'-NaV205 is dominated by an electric dipole active 
mechanism. We have observed the resonant enhance- 
ment of the singlet to triplet transition close to the c- 
axis 68 cm" 1 mode. We described the enhancement of 
the electric dipole transition with the theory of the dy- 
namic DM mechanism and assigned the 68 cm" 1 mode 
to a c-axis optical phonon. From the analysis of the 
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phonon infrared spectra we concluded that two types of 
ladder planes exist with the zigzag charge order along the 
ladders where the charged stripes across the ladders are 
aligned approximately in [110] or in [110] directions. 
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